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Abstract 

The slow passage through a Hopf bifurcation leads to the delayed appearance of 
large amplitude oscillations. We construct a smooth scalar feedback control which 
suppresses the delay and causes the system to follow a stable equilibrium branch. 
This feature can be used to detect in time the loss of stability of an ageing device. As 
a by-product, we obtain results on the slow passage through a bifurcation with double 
zero eigenvalue, described by a singularly perturbed cubic Lienard equation. 
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1 Introduction 



The main motivation for this paper is the problem of ageing of dynamical systems. Assume 
that the behaviour of a device can be described by an n-dimensional ordinary differential 
equation x = f{x). The system is designed in such a way that x* is an asymptotically 
stable equilibrium point, which corresponds to the desired behaviour. Slow changes of 
the system's characteristics due to the ageing process can be modelled by a slowly time- 
dependent equation 

dx 

— = /(x,et), 0<e«l, (1) 

where /(x,0) describes the dynamics of the brand new device. 

Equation (|^) being a nonautonomous differential equation which is difficult to solve, 
one is tempted to consider instead the 1-parameter family of autonomous systems 

dx 

— = /(x,A), A = constant. (2) 
dt 

One hopes that if the "quasistatic approximation" (|2|) has a family of attractors depending 
smoothly on A, then solutions of (0) should be close, at any given time t, to the attractor 
of (I) with \ = £t. 

This is at least partially justified by the following result [PR, VBK| , |F^] : if x*(A) is 



a family of asymptotically stable equilibria of (|2|), then any solution of (||) starting in 
a sufficiently small neighbourhood of x*(0) will, after a short transient, track the curve 
x*{£t) at a distance of order e. For the ageing device, this implies that we need not worry 
as long as the "nominal" equilibrium x*{\) remains asymptotically stable. 

This naturally raises the question of what happens if the equilibrium x*{\) undergoes 
a bifurcation at A = Aq. Such problems are usually referred to as dynamic bifurcations 



I Ben I ■ Bifurcations with a single zero eigenvalue have been studied in some detail. It 



turns out that saddle-node bifurcations will result in a sudden jump of the solution [Hab], 
which may have catastrophic consequences for the device. Transcritical and pitchfork 
bifurcations generically result in a smoother behaviour, where the trajectory follows one 
of the stable equilibria created in the bifurcation ||LS|] . This feature might be used to 
detect the loss of stability of the nominal equilibrium, in order to switch off the device 
before any harm is done. 



The case of a Hopf bifurcation has been analysed more recently |Sh| , Nel, Ne2]. The 
surprising phenomenon is that instead of directly tracking the limit cycle created in the 
bifurcation, the trajectory remains close to the unstable equilibrium for some time, be- 
fore jumping to the periodic orbit. Instead of oscillations with a continuously increasing 
amplitude, one thus observes the sudden appearance, after some delay, of large amplitude 
oscillations (Fig. ||a). The bifurcation delay is stable with respect to smooth deterministic 
perturbations. For the ageing device, this phenomenon implies that the loss of stability 
cannot be detected soon enough to avoid catastrophic oscillations. 

To avoid such problems, one may try to control the system. A simple affine control of 
the quasistatic system (^) would be 

^ = /(x,A)+6n, (3) 

where 6 is a given vector in M" and u is a scalar function. Two cases have been stud- 
ied: 
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1. Open loop control: u{t) is a function of time taking values in some compact interval 
U. The analysis of (|^ for all possible functions of this type leads to the notion of 
control sets, which tend to form around invariant sets. The A-dependence of control 
sets near bifurcation points has been studied in |CK, CHK, 



2. Feedback control: u{x{t)) is a function of the state of the system. The bifurcation is 
"stabilized" by choosing the function u{x) in such a way that the nominal equilibrium 
of (^) is stable when A = Aq. For topological reasons, this makes the bifurcation 
supercritical, and should avoid exploding trajectories at least for A slightly larger than 
Ao [Fi, N,|MS|. 



The sometimes surprising behaviour of dynamic bifurcations shows that we should 
consider the slowly time-dependent version of (^) as well. Here we will limit ourselves to 
feedback controlled systems of the form 

^ = f{x,et) +buix,et) (4) 
at 

in the specific case where / undergoes a Hopf bifurcation. Since we wish to analyse 
on the time scale e~^, we introduce the slow time t = et and rewrite @) as the singular 
perturbation problem 

dx 

= fi^^T) + bu{x,T). (5) 

Our aim is to design a feedback control u in such a way that the bifurcation delay, and 
hence the sudden appearance of large amplitude oscillations, are suppressed. 

A feedback affecting only nonlinear terms m x — x* will have no effect on the delay. 
We thus have to modify the linearization ^ of / at the bifurcation point. Shifting the real 
part of the eigenvalues of A will merely postpone the problem to some later time. The 
only solution is thus to shift the imaginary part of the eigenvalues in order to produce a 
double zero eigenvalue. 

Our strategy will thus be the following. First we design a feedback in such a way that 
when the parameter A is varied through the bifurcation value Aq, new stable equilibrium 
branches are created at a distance of order (A — Aq)^''^ of the nominal branch x*(A). The 
resulting vector field turns out to be a codimension four unfolding studied in |lKKR[| . We 
then show that the solutions of the corresponding system (|5|) actually track one of these 
branches, a feature which can be used to detect the bifurcation point. 

To analyse the time-dependent system, we extend certain methods of [|B1| (see p2(| 
for a summary and |BK] for applications). However, this paper is written in a largely 



self-contained way. It is organized as follows: the main result is stated in Section ^ in 
Section ^ we recall Neishtadt's result on bifurcation delay | |Nel ], which will guide us in 



the construction of the feedback in Section ^. Section ^ is devoted to the proof of the 
main result for the singularly perturbed equation (|5|) . Finally, in Section ^, we give a few 
remarks on what happens when the control is slightly imperfect, and the requirements of 
the theorem are no longer met. 

Acknowledgments: It is a pleasure to thank Klaus Schneider for welcoming me at the 
Weierstrafi Institute and for inspiring discussions on this interesting topic. I thank Dima 
Turaev for pointing out useful references on high codimension bifurcations. This work is 
supported by the Nonlinear Control Network of the European Community, Grant ERB 
FMRXCT-970137. 
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2 The Problem and Main Result 



We consider the feedback controlled dynamical system 
d T 

— = f(x,\) +bu(x,X), xeM", AeM (6) 
at 

where: 

• The uncontrolled vector field f{x,X) undergoes a Poincare-Andronov-Hopf bifurcation 
at the origin. 

• The vector 6 G M" is imposed and fixed. It describes the direction in which the system 
can be steered. 

• The scalar feedback control u{x, A) G M is a function to be determined in such a way 
that the solution behaves "smoothly" when A is slowly varied. Its dependence on x 
and A should be as simple as possible (e.g. polynomial). 

More precisely, we will assume that the uncontrolled vector field satisfies the following 
hypotheses: 

(HI) Domain and smoothness: The function f{x,X) : P x / ^ M" is analytic in a 
neighbourhood T> of the origin in R" and an interval / containing 0. 

(H2) Hopf bifurcation: There exists a curve a;* (A) : / — > M" with x*(0) = and 
/(x*(A),A) = 0. The Jacobian matrix A(X) = d.j:f(x*{X),X) G M"^" admits two 
eigenvalues a(A) ± ic(;(A), where a(0) = 0, a'(0) > and a;(0) = / 0. All other 
eigenvalues of ^(A) have a strictly negative real part. 

These hypotheses imply in particular that there is a set of coordinates x = {y, z), with 
y G M"" (m = n - 2) and z = (^,77) G such that 

fix, 0) = (^-y + 'A , Ao=( ^ , (7) 

where all eigenvalues of A- g have negative real part and and go are of second 

order in y and z. In these coordinates, we write b = ( j,~ ) ■ 

(H3) Controllability: bo / 0. 

Since Aq is rotation invariant, we may assume that 60 = (i)- 

The next hypothesis is more technical and its meaning will become clear in Section ^. 
It is, however, generically satisfied. 

(H4) Nondegeneracy: Define the matrix T = —{AZ^b-,Az'^b^) G IR*"^^. Let g^ denote 
the first component of go, and h{z) = g^{Tz,z). Then either d^^h{0) 7^ or 
%^/i(0) < 0. 

Finally, we require the following property of the equilibrium branch x*(A): 
(H5) Velocity of equilibrium: The //-component of dAX*(A)|^^Q is different from zero. 
The main result of this paper is the following. 
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Figure 1. (a) The slow passage through a Hopf bifurcation (occurring at t = 0) leads to 
the delayed appearance of oscillations, (b) The control we construct suppresses this delay, 
and causes the orbit to track a curve lying at a distance of order r^^^ from the nominal 
equilibrium (represented by the axis ^ = 0). These plots actually show a solution ^(r) 
of equation ( pi] ) with a(r) = r, j{t) = S{t) = 0, R{^,r],T,s) = 1, e = 0.003 and (a) 
Ii{t) ~ — 0.2, (b) ii{t) — 2.5t. The eigenvalues of the linearization around the origin are 



Theorem 2.1. Assume that the hypotheses (H1)-(H5) hold. There exist 

• strictly positive constants T, M and k, 

• a neighbourhood J\f <ZT) of the origin in M", 

• a smooth feedback control u{x, A) : x [— T, T] — > M with u{x*{X), A) = 0, 

• and a curve x+ (A) : [0,r] ^M" with limA^o+||a;+(A)-x*(A)||/\/A = i^ 

with the following property. For every tq E [— T, 0), there exist strictly positive constants 
c\ and eo such that, if < e < eq, any solution of the equation 



= /(^'^) + bu{x,T) 



(8) 



with initial condition x{to) G Af exists on the interval [to,T] and satisfies the following 
bounds: 



\x{T)-X*iT)\\ ^ M^, 

\t\ 

\x{t)-x*{t)\\ ^ Me^/\ 

-/4 



(r)-x+(r)|KM - + - 



To + cie|lne| ^ r ^ - e^/^, 



(9a) 
(9b) 
(9c) 



This theorem shows that solutions of (^) will first track the nominal equilibrium curve 
x*(t) for negative r, and then track a new equilibrium curve x+(r), situated at a distance 
of order ^fr from x*(r), for positive r (Fig. |l|b). This result looks similar to the stability 
exchange for pitchfork bifurcations in [|LS| . It is, however, more difficult to obtain because 
the specific nature of the bifurcation with double zero eigenvalue makes the problem in- 
trinsically two-dimensional. In particular, solutions tend to rotate around the equilibrium 
branches which are foci near the bifurcation point. 

The proof is divided into two main steps. In the first one, described in Section |^, we 
consider the autonomous system (^. We construct a feedback u(x, A) in such a way that 
after some changes of variables, including a center manifold reduction and a transformation 
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to normal form, the dynamics is governed by the two-dimensional effective equation 
d^=^ 



(10) 



where /i(0) = 7(0) = (5(0) = and 0(4) = 0{\\{^,r])\\'^). The function ^(A) can be 
controlled to some extent by the choice of u. This so-called cubic Lienard equation happens 
to be a codimension-4 unfolding of the vector field (r/, — — C^??) which has been studied in 
detail, see [KKR| and references therein (the authors in |KKR| actually use an equivalent 
unfolding obtained by the transformation C + ^7)- 

The second step is to show that with this particular feedback, the slowly time-dependent 
system (H) can be reduced to 



d^ 
"d7 = ^ 

= /"(^)? + Mr)r] + l{T)e + S{t)^V + 0(4) + eR{t rj, r, e). 



(11) 



The function i?(0,0, r, 0) is related to the drift dx*(r)/dr of the equilibrium branch, and 
Hypothesis (H5) imphes that R{0, 0, 0, 0) / 0. In Section | we prove the following result. 

Theorem 2.2. Assume that /u'(0) > 0. There exist positive constants d, T, M, k and a 
neighbourhood J\A of the origin in with the following property. For every tq G [— T, 0), 
there is a constant ci > such that for sufficiently small e, any solution of (|Tl| ) with 
initial condition iS,,'n){To) ^ satisfies 



\ar)\^M- 



\rj{T)\ ^ M- 



|r|i/2' 

|e(r)| ^ Me^^ |77(r)| ^ Me^/s, M -(|) 

where Ti(e) = tq -|-cie|lne|. If, moreover, the relations 

a'(0) 



for Ti (e) ^ T ^ 
e\2/3 



£\2/3 
d) ' 



, e\2/3 



1 



i?(0, 0,0,0) / 



(12) 
(13) 

(14) 



hold, then for (e/d)^/^ ^ r ^ T we have 

\vir)\<M 



S £ 
- + 



1/2 



r rV4 

^ +£l/2^1/4g-«rV. 



r 



1/2 



(15) 



where 



e+(r) 



^ + 0(r), i/i?(0,0,0,0) >0, 
-^ + 0{t), i/i?(0,0,0,0) <0 



(16) 



are equilibria of (|To|), i.e., the right-hand side of (|T^) vanishes when ^ = ^_|_ and r] = 0. 
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3 The Uncontrolled Problem 



We state a version of Neishtadt's result on bifurcation delay for the uncontrolled problem 

e^ = /(x,r). (17) 



The results in |Nel, Ne2| apply in fact to a more general slow-fast system. We are 
interested in the computation of the delay time. To this end, we define the function 



^{t)= p(s)ds, p{s) = a{s) + iu;{s). (18) 
Jo 

In some interval (— T, 0) in which a(ro) < we can define the map 

n(ro) = sup{r I Reip{s) < ReV'(ro),To < s <t}, (19) 

r>ro 

giving the first time at which ReV'(r) becomes equal to Re'i/'(7"o) again. Observe that 
n(ro) > 0, lim^„^on(r) = and lim^„^o n'(ro) = -1. 

Theorem 3.1. Assume that /(x,A) satisfies Hypotheses (HI) and (H2) for all (x, A) in 
some complex neighbourhood ofD x I. Let tq < be such that a(r) < for tq ^ r < 0. 
There exist a neighbourhood N of x*{tq), a constant M > 0, a buffer time r_|. > and a 
continuous function (5(e) going to as e — > 0, such that any solution of ( p^ with initial 
condition x{tq) £ J\f satisfies 

\\x{T)-x*{T)\\i^Me for To + 5{e) i^T s^f -6{e), (20) 

where the bifurcation delay time f is given by 

f = min{n(ro),r+}. (21) 

The quantity f gives in fact a lower bound on the bifurcation delay. Under more strict 
assumptions, it also gives an upper bound in the limit e — > 0, see [ Ne2(| and the article by 
F. and M. Diener in [Ben]. 

The important fact for us is that f is the minimum of two quantities. The first one, 
n(ro), depends on the initial condition and thus cannot be used to influence the delay. 
We thus have to modify the buffer time r+ , which is determined in the following way. The 
function ^'(r) and the solutions of ( p^ can be continued to a complex neighbourhood of 
r = 0. For sufficiently small |ro|, the real times tq and n(ro) can be connected by a path 
r(ro) (lying in the upper half plane if ujq < 0), along which Re^'(r) is constant. One 
defines a negative buffer time t_ as the smallest real time such that r(r_) exists and has 
certain properties stated in |Ne2|. The positive buffer time is given by 

T+ = sup n(r). (22) 

T_<T<0 

The existence of this buffer time is a rather subtle, nonperturbative effect, which cannot 
be understood by naive perturbation theory. The proof uses an integration of ( p!7[) along 
a path r(r). 

Example 3.2. If a(r) = r and uj{t) = —loqi, then the level lines of 

Re^(r) = ^[{Rerf - {luiT - Loof + io^] (23) 

are hyperbolas centered at r = iujq. Thus we have n(ro) = —tq and the buffer times are 
T± = iojQ. The delay may be suppressed by choosing a control in such a way that ujq = 0. 
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4 Design of the Feedback Control 



We start by analysing for A = 0. Coordinates are chosen as in where we can scale 
time in such a way that ujq = 1. Our feedback control will be constructed to satisfy two 
requirements: 

1. The properties of the bifurcation delay time show that we should move the imaginary 
part of the eigenvalues ibu;(0) towards the origin, in order to shift the buffer time to 
the bifurcation time. 

2. This will produce a bifurcation with double zero eigenvalue. In analogy with works 



on stabilization of bifurcations |Ae, MS], the equilibrium should be stable at the bi- 
furcation point, in order to avoid escaping trajectories. 

We claim that an appropriate feedback control has the form 

u{x,0) = {1- i^)S,+v{z). (24) 

[We recall that x = {y, z) G M™' x and z = (^, ry).] The parameter u will be ultimately 
set to 0, we introduce it in order to keep track of the effect of an imperfect control. Its 
aim is to shift the eigenvalues of the linearization to the origin. The function v{z) 

is a nonlinear term which should assure that the origin of (^) is stable. 

4.1 Center Manifold Reduction at A = 

With the feedback control (pi), equation (pD takes the form 



y = A-y + g^{y,z) + B-Z 
z = Aqz + gQ{y,z), 



(25) 



where g- = g-{x) + bv{z), go = go{x) + bv{z), = (1 - i^){b-,0) e M™^^ ^j^g 
marginally stable part of the linearization is given by the matrix Aq(i') 



1^ 

-u 0, 



Proposition 4.1. Equation ( p5| ) admits an invariant manifold on which the dynamics is 
governed by the equation 

z = Aoz + G{z), (26) 



where G{z) = 0{\\z\\ ) is given to third order by equation (31) below. 



Proof: The existence of the manifold follows quite directly from the Center Manifold 
Theorem |Ca|] . We would now like to compute G. First, we introduce the matrix T G M'"^^ 
satisfying A^T — TAq = —B^, which exists because A_ and Aq have no eigenvalues in 



common [Kr, Wa|. In fact, it is given by 

T = {ti,t2), ti = A_t2, t2 = {u - l){A'i + uty^b-. (27) 

The change of variables y = yi + Tz yields 

m = A^yi +9^{yi,z), g-{yi,z) = g^{yi +Tz,z) - Tgo{yi +Tz,z), 

z = Aqz + go{yi,z), go{yi, z) = go{yi + Tz, z). (28) 
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This system admits a center manifold locally described by yi = h(z), where h{z) satisfies 
the partial differential equation 

A-h{z) + ~g-{h{z), z) = dMz) [Aoz + go{h{z), z)] . (29) 

For a vector field F{x) : M" M™, we denote by F^''\x) the terms of order k of its 
Taylor expansion around 0. We know [Ca| that h{z) = h^'^\z) + 0{\\z\\^), where 



d,h^'^\z)AQZ - A^h^^\z) = ~g^^\Q,z). (30) 
The motion on the center manifold is thus given by (p6|), where 



G{z)=go{h{z)+Tz,z) 

= G^^\z)+G^^\z) + 0{\\z\\^), 
G^^\z)=g^^\Tz,z) ^''^ 
(z) = 5^ {Tz, z) + dyg'i^ {Tz, z)h(''^ (z). 



To compute G^^\z), we need to solve equation ( |30|) for h^'^\ In fact, we will only need to 
know that h^'^\z) = —Az}g^i^ + 0(?y), where g^ is the coefficient of in ^^^^(0, z). □ 

4.2 Stability at A = 

We intend to construct the nonlinear part of the feedback control of the form 

v{z) = vii'^ + V2ilt] + W3??^ + Vif, (32) 

with appropriate coefficients Vi. They will be determined by the following result of normal 
form theory: 

Lemma 4.2. Consider the system 

i = r] + ciC^ + c2Cv + C3f?^ + Cif + 0{r]\\zf, \\zf) 
= die + d2ir, + d^if + ^4^3 + 0{r^\\zf, ). 

Let a = cids + C4 and 13 = 2c^ + d^. Then 

• If di = d2 + 2ci = and a, /3 < 0, the origin is asymptotically stable. 

• Conversely, if the origin is stable, then di = d2 + 2ci = and a, P ^ 0. 



(34) 



(33) 



Proof: The normal form of the 3-jet of ( |33| ) can be written as 

i = v 



with 'J = d2 + 2ci and 6 = di. The assertion has been proved in |Ta|, see also |GH]. In 
fact, the origin is an unstable Bogdanov-Takens point if 7 7^ or if (5 7^ 0. If 7 = 5 = 0, 
it is asymptotically stable if a and /? are both negative, and unstable if one of them is 
positive. □ 

We would like to choose the coefficients of v{z) in such a way that the vector field (|26 



satisfies Lemma 4.2. The next lemma shows that this is generically possible. 
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Lemma 4.3. Assume that the component of go{Tz, z) has the /orm 01^^+04^^+0(7/, ^^). 
Then, if either ci 7^ or C4 < 0, one can find a function v{z) of the form (^) such that 
the origin in (^) is asymptotically stable. 

Proof: Let us consider 

The first two conditions of Lemma |4.2| are satisfied if we clioose vi = —di and V2 = 
— (^2 + 2ci). To satisfy the third one, we have to choose ^3 in such a way that a = 
ci(c?3 + V3) + C4 < 0, which is possible under our assumptions. The last condition looks 
more difficult to satisfy, but in fact, we have f3 = U4+constant, where the constant depends 
only on previously fixed quantities, so that /5 can always be made negative. □ 

The requirement on gQ{Tz, z) is nothing but Hypothesis (H4). Note that the coefficients 
U3 and ^4 only have to satisfy inequalities, while vi and V2 must have a specific value. We 
will discuss in Section ^ what happens when these coefficients are not exactly equal to the 
prescribed value. 

4.3 Choice of u{x, A) 

We consider now equation (P) for general, fixed values of A. If x*(A) is the equilibrium 
branch of /, an affine transformation x = x* (A) + S{\)[y) yields the system 

y = A_ {\)y + 5_ (y, z, A) + 6_ {\)u{y, z, A) ^^^^ 
i = Aq{X)z + 5(o(y, z, A) + bo{X)u{y, z, A), 

where ^-(A) has eigenvalues with negative real part for sufficiently small A, and 

In fact, since ^0 is rotation invariant, we may assume that 60(A) = (5). The terms g- 
and go are nonlinear. 

We now choose the feedback control u{x, X) in such a way that 

u{y, z,X) = {l + C A) [(1 - u)^ + v{z)] , (38) 

where C is some constant to be determined, and v{z) has been constructed in the previous 
section. In this way, the linear part of i becomes 

^'^^^ = {-u;{X) + {i%{l + CX) Iw)- ^^^^ 
This matrix can be further simplified by a shearing transformation: 

where 

;,(A) = -a;2 + (1 _ + C7 x)uj = (C - uj'{0))X -u + 0{X^ + u^). (41) 
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Figure 2. Schematic bifurcation diagram of equation (43a) in the plane 7 = (5 = 0. The 
hne A-B is the original Hopf bifurcation. By moving the eigenvalues' imaginary parts to 0, 
we change the function /i(A) in such a way that /i(0) = 0. This produces new bifurcation 
lines. The transition A-F is a supercritical saddle-node bifurcation, the transition C-B a 
subcritical one. D-C is a subcritical Hopf bifurcation, D-E a homoclinic bifurcation and 
E-F a saddle-node bifurcation of periodic orbits. 



In fact, any of the pairs (A, u), (A, fi) or (a, /x) can be considered as independent parameters 
used to produce the bifurcation with double zero eigenvalue. We will henceforth set u = 
and consider (a(A),;u(A)) = (C — u;'(0))A -|- 0(A^)) as a path going through the origin of 
the two-dimensional parameter space. 

The system ( pG] ) admits a center manifold described locally by its parametric equation 
z = h{y, A). On this manifold, the dynamics is governed by the equation 

z = B{X)z + G{z,X). (42) 

It will not be necessary to compute G{z, A). We only need to know that G{z, 0) has been 



computed in Proposition 4.1 and satisfies the requirements of Lemma |4.2|. 



4.4 Normal Forms 

The nonlinear term G{z, A) can be simplified by eliminating all terms which are not 
resonant at A = 0. A convenient basis of resonant terms up to order 3 is given by 



Proposition 4.4. The 3-jet of (^2|) can be written in either of the following forms: 

■q = IJii + at] + + 5^7] - ^ - i rj, 
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^ = 7] 



(43b) 



(43c) 



Each set of parameters (a, ^, 7, 5), {a, (3, 7, 6) or {fJ-i, f^s, ^^a) depends on A and vanishes 
at X = 0. 



Proof: (43a) is obtained by eliminating nonresonant terms and rescaling space and time. 
( |43b| ) comes from the linear transformation 7/ 1— > 77 — followed by a scaling. ( [43cj ) is 
obtained with the transformation ^ ^ ^ + ^7. □ 

These normal forms describe equivalent codimension-four unfoldings of the singular 
vector field (^ = r],rj = —frj — f). Indeed, the bifurcation is produced by setting four 
parameters to zero: the initial bifurcation parameter A, the imaginary eigenvalue shift 
and the coefficients vi and V2 of the feedback control. 

Equation ([43bD has been studied by Takens in the particular case 7 = (5 = (which 
occurs when there is a symmetry z 1— > —z), who obtained the bifurcation diagram shown 
in Fig. § [Ta, GH |. A codimension-three unfolding corresponding to /i4 = has been 
studied in VT |. The full codimension-four unfolding (43c) is discussed in | KKR[| (see also 
references therein). 



5 Slowly Drifting Parameter 

Let us now consider the slowly time-dependent equation (P). We write /(x, \) + bu{x, A) = 
F(x,X)^ where u{x,X) is the feedback control constructed in the previous section. The 
equation thus reads 

e— = F{x,T). (44) 

By construction, F{x,t) vanishes on the original equilibrium branch x*{t). Let us write 
F{x*{t) + xi,t) = A{t)xi + g{xi,T), with g = 0{\\xip). Hypothesis (H2) and the 
properties of u imply that there exists an interval [— Ti,0) in which the matrix A{t) has 
only eigenvalues with a negative real part. Moreover, these real parts are bounded away 



from on any interval [ro,ri] C [— Ti,0). It is well known (see for instance |PR|) that 
this implies that any trajectory starting at tq in a sufficiently small neighbourhood of the 
origin reaches an 0(e)-neighbourhood after a time of order e|lne|, where it remains until 
r = Ti (this result is proved using Lyapunov functions) . Hence it suffices to study ( ^4|) for 
r ^ Ti, with an initial condition x{ti) = 0(e). 

We first simplify equation (^^ by applying similar transformations as for the time- 
independent system, which will produce some additional terms of order e. We will obtain 
an effective two-dimensional equation, which we will study by various methods in three 
different time intervals. 
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5.1 Reduction of the Equation 

The time-dependent translation x = x*{t) + xi yields the equation 



dxi 



A{t)xi +g{xi,T) - £ 



dx*{T) 



(45) 



dr dr 

For small r, there exists a nonsingular matrix S{t) such that S^'^AS is block-diagonal, 
with blocks A_ and B as in ( p5| ) and (42). The transformation xi = S{t)x2 yields the 
equation 



dX2 

' d7 









B{t) 



X2 + S g{Sx2,T)-eS 



-1 



d5(r) ^ ^ dx^(T) 
dr dr 



(46) 



Next we apply a center manifold reduction. We write X2 = {y2, Z'l) G I^*" x 1^^- 

Proposition 5.1. For sufficiently small e , equation (|4^ ) admits a local invariant manifold 
(local in Z2 and t), with a parametric equation of the form X2 = {h{z2,T,£), Z2), where 
11^(2:2, r, e)|| ^ M(||z2|p + t'^ + £)■ Any solution starting at a distance of order e from this 
manifold is such that y2{T) = h{z2{T),T, e)+0{e) in some neighbourhood of {x2,t) = (0, 0), 
where Z2{t) satisfies the equation 

dz2 



dr 



with B{t) and G{z2,t) the same functions as in 
of P{0, 0, 0) is different from 0. 



B{t)z2 + G{z2,t)+£P{z2,T,£), (47) 

1). Moreover, the second component 



can also be written as 

^_(r)y2 + ew^{T) + 0(||y2f + 1^2^ + r^ + e^) 



Proof: Equation (^ 

<^ty2 

dtZ2 = B{T)z2+£WQ{T)+0{\\y2f + \\z2f+T'^+e^ 
dtr = £ 
dte = 0. 



(48) 



The transformation y2 = ys — eA_^(0)w_(0) yields a system whose linearization at the 
point (y3, ^2, r, e) = (0, 0, 0, 0) is the matrix 



/A_(0) \ 

B{0) wo{0) 

1 

V / 



(49) 



The center manifold theorem implies the existence of a local invariant manifold 1/3 = 
h{z2,T,£) = 0(||z2|P + t"^ +e'^). Moreover, it is shown in Lemma 1, p. 20 of [Ca|, that in 
some neighbourhood of the origin, any solution satisfies a bound of the form 

112/3 (r) - hiz2iT) , r, e) K M e-'^(^-^«)/^ ||y3(ro) - /i(z2 (ro) , r, e) || , (50) 

for some positive M,k. Since ||?/3(ro) — /i(z2(ro), r, e)|| = 0{e), we have 

= B{t)z2 + goih{z2,T,£) + ©(e), Z2, r, e). (51) 

dr 

When £ = 0, this equation coincides with the equation on the instantaneous center mani- 
fold. The assertion on P(0, 0,0) follows from Hypothesis (H5). □ 
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As a final reduction step, we may apply a nonlinear transformation putting the 3-jet 
of G{z2,t) into one of the canonical forms (43). This will produce terms of order e||2:2|p, 



that we may also absorb into the remainder P{z2,t, e). A redefinition of the variable rj 
yields the canonical form In the sequel, we will use the equivalent form 



(52) 

£7) = /.(r)C + a(r)7? + 7(r)^2 ^ _ ^3 _ ^2^ ^ 0(||^||4) ^ ^^(^^ 

where we now use the dots to indicate the derivative with respect to r. We can rescale 
time in such a way that 

a{T) = cT + OiT^), Mt) = t + 0(t2), j{T) = Oir), 5{t)=0{t), (53) 

where c = da/d^(0) = d(0)/(C - u;(0)), see 

We will study equation (^) by different methods in three different regions. They are 
characterized by a constant d which will be chosen sufficiently small, but is independent 
of £. These regions are: 

1. Before the bifurcation: ri(e) ^ r ^ — (e/d)^/^, where ri(e) = tq + O(elne) is chosen 



in such a way that z{ti) = 0{e). In this region (52) can be transformed into a one- 
dimensional complex equation, that we study by suitable normal form transformations. 
The main difficulty is to use the averaging effect of fast oscillations around the focus. 

2. During the bifurcation: — (e/d)^/^ ^ t ^ (e/d)^/^. In this inner region, ( |5^ ) can be 
reduced to a time-dependent Hamiltonian system by an appropriate scaling. 

3. After the bifurcation: (e/d)^/^ ^ r ^ T, where T will be chosen sufficiently small. By 
choosing c in an appropriate way, we ensure the existence of an attracting equilibrium 
branch. 

5.2 Before the Bifurcation 

In the first region, we use a coordinate transformation which diagonalizes the linear part 
approximately. It is given by 

C = [e-/Vl'/'C + , (54) 



and its inverse reads 



(55) 



V2l 

In these variables, equation ( ^2|) becomes 

eC = p{r)C + e^{T)C + ew{T) + Go(C, C, r, e), (56) 

where 

p{t) = a{T) + i + eV(r), V'(r) = 0{\r\-'/^), 

(^(T)=iA(T)/4/.(r)=0(|T|-i), 



w{t) = -j= e- i^/Vr '/'^(0, 0, r, 0) = 0(|t|-i/^), 

Go(C,C,r,e)=0(|T|i/^|C|2) + 0(|r|-^|C|3) + 0(|r|-5/4|Cn. 
The main result of this subsection is the following estimate. 



(57) 
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Proposition 5.2. // d and e are sufficiently small, there exists a constant Mi > such 
that equation ( |56[) admits a particular solution Co{t) satisfying 



Co(r)+e- 



w{t) 



^ Ml 



-|9/4 



forn^^Ti^ - (58) 



a(r) + i y^-/i(r) 

Moreover, any solution of ( |56D wzi/i initial condition C(''"i) = ^(s) satisfies the bound 

\C{r)-Co{r)\^M2ee-<-"-^?^/'' (59) 
on i/ie same izme interval, for some positive constants M2, n. 



Corollary 5.3. Any solution of (^) with initial condition of order e satisfies 

e(T)=0(e|T|-i), ii{r) = 0{e\T\~^l^) 
on the interval ti ^ r ^ — (e/d)^/^. In particular, at t = — (e/d)^/^ we have 
^ = d2/3ei/3 [^(0, 0, 0, 0) + 0{d) + 0{{e/d)^/'^)] , 



(60) 



(61) 



Remark 5.4. One can in fact show the existence of a particular solution admitting an 
asymptotic series of the form 



Co(r) 



-|3/4 L 



Co{t) + Ci{t)j-^ + C2(t) 1^ + • 



(62) 



For r = (e/d)^/^, we get an asymptotic series in d. Proposition 5.2 is only a first step, but 
suffices for our purposes. 



The proof of Proposition 5.2 is based on the following two lemmas. The first one is a 



rather trivial, but very useful bound on an integral we will encounter several times, while 
the second one gives an "a priori" estimate on Coi'^)- 

Lemma 5.5. Assume that we are given 

• constants ri ^ r < 0, oq, wq > and p, q such that p + 1 — q > 0. 

• ^'(t) : [ti,0) C differentiahle such that Re^'(r) < and |^(r)| ^ \T\P/aQ. 

• w{t) : [ti,0) C differentiahle with \w{t)\ ^ wo\t\'^~^ and \w{t)\ ^ wo\t\'^~'^. 
Then 



C 

^ 1111 — K, K = qqWo 



2 + 



1 



p+l-g 



Proof: Use integration by parts once. 



(63) 
□ 



Lemma 5.6. Let cq > be given. If e and d are small enough, there exists a constant 



Mi(co) such that any solution of (56) with initial condition |C(''"i)| ^ cqS satisfies 



IC(t)| ^Mi- 



-|3/4 



forn^Ti^ - (e/(i)2/3. 



(64) 
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Proof: Notice that for small enough d, we have \p\ ^ i^ol''"!^^^- 

• Step 1: Simplification of the linear part. 
Consider the initial value problem 

ss = [p{t) - p(t)] s + eifir) - elp{T)s^, s{ti) = 0. (65) 

Since drlsp = {ips — Tps){l — |sp), its solution satisfies \s{t)\ ^ 1. (In fact, one can 
prove that s = 0{£\t\~^^'^)). The transformation C = Ci + siT)Ci yields 

eCi = Pi(t)Ci + ewi(r) + Gi(Ci,Ci,^,e), (66) 

where pi = p + eips satisfies \pi\ ^ Kilrl-*^/^, and Gi and wi satisfy similar bounds as 
Go and w. 

• Step 2: Simplification of the cubic part. 

Let us consider the effect of a change of variables 

Ci = C2 + MC2,C2,t), MC2,C2,t)= Y1 ^nm{r)Q'^. (67) 

n+m=3 

It transforms ( |6^ ) into 

(2) 

^C2 = P1C2 + — edc^hwi — ed-^hwi + + 

+ [Gf + pih - d^hpiC2 - 5^/ipiC2 - ^drh] + 0{\\z2f), (68) 

(k) 

where G\ denotes terms of order k. In particular, 

G?{C2,C2,r)= Yl 9nmir)C^C7, (69) 

n+m=3 

with gnmi'T) = C(|''"|~^)- We see that the term in brackets of (|68| ) can be eliminated 
if 

eKm = [(1 - n)pi - napi] hnm + 9nm- (70) 

This is a linear equation which can easily be solved. If (n, m) ^ (2, 1), one can choose 
the initial condition in such a way that \hnm{T)\ ^ K2\t\-^/'^. When (n,m) = (2,1), 
however, we cannot obtain such a good bound because the imaginary part of the term 
in brackets vanishes. Thus we do not attempt to eliminate this term. We obtain the 
equation 

eC2= [pi{t) - g2i{T)\C2\^]C2 + ewiir) + G2{C2X2,r,e), (71) 

where 

|G2(C2,C2,r,.)| ^M2[|r|VY2|2+e|^|-7/4|^^|2 + |^|-5/4|^^|4]_ (73) 

Let us write />2(t, C2) = Piir) - 5'2i(T")|C2p- 
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• Step 3: Proof of the hound (64). 
The solution of ( |7l[ ) has to satisfy 

C2(r)=e[*^«-*^(-^)]/^C2(ro) + 

e*2(r)/. re-i'2M/.[^^(^) + lG2(C2,C2,^,e)] ds, (73) 

where ^2{t) = Jq P2{s, C2{s)) ds. We define the time 

T* = sup {r I |C2(s)| ^ for n ^ s ^ r}. (74) 

[ri,-(e/d)2/3] 

By continuity, r* > ri. Using the bounds on G2 we can show that for ri ^ r ^ r*, 

\P2{tX2)\^K2\t\^/^ 

\G2\^Ms[e\T\'/' + eW-'/'], (75) 
|d.G2KM3[e|rr^/^ + e2|rrii/^], 

where dT-G2 = drG2 + ^2C2 + d-^JJ2C2 is estimated using ([7l|). Using Lemma 5.5 
and the bounds on wi, we obtain that 

|C2(r)KM4[.|r|-3/4+eV|-9/^] 
^ |C2(r'^)KM4ei/2[di/2 + d3/2]. ^ ^ 

Taking d sufficiently small, we have IC2 (''"*) I < £1/2. If we assume that < -{e/df/-\ 
we contradict the definition of r*, which shows that r* = — (e/d)^/^. 

Going back to the initial variables and using the bounds on s and /i, we obtain the 
conclusion of the lemma. □ 



Proof of Proposition |5.2| . After subtracting £w{t)/p{t) from ( |56D and eliminating 
the term linear in C, in the same way as in Lemma |5.6|, we obtain 



sQi = pi{t)Ci + e^wi{T) + Gi{CiXi,r,e), (77) 

where 

\Pi{t)\ ^Ki|t|1/2, 

wi{t)=0{\t\-'I^), (78) 
Gi(Ci,Ci,r,e)=0(|T|VYi|2) + 0(e|Tr^/4|CiP) + 0(|T|-i|Ci|^). 



Let ^i(t) be a primitive of Pi(t). The solution of ( [77| ) with initial condition Ci(''"i) = 
should satisfy 

Ci(r) = e*iW/^ re-*i(^)/^[eu;i(s) + 1gi(Ci, Ci, e)] ds- (79) 



Using the a priori estimate of Lemma |5.q and the bounds (|78D, we can apply Lemma 5.5 
to estimate the integral and obtain (|58|). To obtain (|59|), it is sufficient to subtract the 
particular solution Co(''") from the general solution, and to use the modulus as a Lyapunov 
function. □ 
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5.3 During the Bifurcation 

In this subsection, we study (^) on the time interval [— (e/d)^/^, (e/d)^/^]. We recall 
that d is a constant which will be chosen small, but is independent of e, while e is small 
with respect to d. In other words, we consider the system on a time scale d~^, which is 
intermediate between 1 and e~^. In fact it turns out to be useful to take e ^ d^. 

Prom Corollary 5.3, we know that at the time r = — (e/d)^/^, ^ is of order e^/^ and r/ is 



of order e^/^. The basic idea to analyse the motion during the bifurcation is to introduce 
the scaling of variables 

e = ie/d)'/'x, rj = (e/d)2/3y, r = {e/df^H. (80) 

This scaling transforms the system (^) into 

di: = y + {s/d)^/^Qi{x,y,t,e), 

dy = tx-x-^ + dRo + {e/dY''^Q2{x, y, t, e), 

where Rq = R{0, 0, 0, 0) and the functions Qi and Q2 are uniformly bounded by constants 
independent of e and d. This system should be studied on the time interval t £ [—1,1], 
and with an initial condition 

x(-l) =d[Ro + 0{d) + 0{{e/d)'/^)] = dRo + 0{d^), 
y{-l) = d [0{d) + 0((£/d)i/3)] = 0(d2). 

Equation (^) is a small perturbation of the Hamiltonian system 



1 



d 



1 2 1^4 1^^2 



H{x, y,t) = - -y^ + -x'^ - -tx^ - dRox . (83) 



2" 4 



Lemma 5.7. For — 1 ^ t ^ 1, the solution of equation ( |8l| ) and of the Hamiltonian 
system ( |83|) with the same initial condition differ by a term of order (e/d)^/^. 



Proof: One can use a standard averaging result, see for instance |GH|| , Theorem 4.1.1. 
page 168. □ 

The Hamiltonian (p^) describes the motion of a particle in a potential V{x, t) = 4x'^ — 



4" 

^tx"^ — dRox which changes from a single well to a double well as t grows from —1 to 1 



2' 

(Pig. I). Note that (H) imphes that the energy is negative at t = — 1 if d is sufficiently 
small. Since 

-^{x, y, t) = — (x, y, t) = -^<tf, (84) 

the energy must decrease. Since the barrier between the potential wells always has positive 
energy, the particle has to remain in the same well, which is the right one if i?o > and 
the left one if i?o < 0. 

We will however need a more precise estimate, showing that the orbit remains close to 
the bottom of the well. This is the main result of this section. 

Proposition 5.8. Assume that Rq 7^ 0. For sufficiently small d, the solution of the 
Hamiltonian system ( |83|) with the initial condition (^) is such that 

x{l) = signiZo + 0{d), y{l) = 0{d). (85) 
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Figure 3. The dynamics near the bifurcation point can be reduced to the motion of a 
particle in a time-dependent potential which transforms from a single to a double well. 
The potential V{x,t) = jx"^ — ^tx^ — dR^x is shown for dRo = 0.1 and t = —1, and 
1. Since the energy is decreasing, and the barrier between the wells always has a positive 
energy, the particle remains in the same well. 



Together with Lemma this implies 

Corollary 5.9. During the time interval [— (e/d)^/^, (e/d)^/'^], the solution of ( |52D with 
initial condition given by Corollary \5. satisfies ^ = 0{e^/^) and i] = ©(e^/^). At the time 
T = (e/d)^/^ we have 

^ = (£/d)i/3 signing + 0(ei/3rf2/3) ^ o{{e/df^^), 

ri = 0{e^/^d^/^) + 0{e/d). ^^^^ 



Proof of Proposition |5 . 8| . To simplify the notation, we consider the case Rq = 1. 
• Step 1: Transformation of the equation. 

Let Xi,{t) be the positive solution of tx^ — + c? = 0. Then we have 

Mt) = xS) = (87) 

a{t) = 3x1 -t^ max{-t, Id"^/^, 2t}. (88) 
We decrease the order of the drift term d by the transformation 

X = xi + Xi,{t) — d^ — 
y = yi + dx^(t) - d^-^[-^)^ 
which implies that 1), 1) = 0{d'^) and yields the system 
dxi = yi 

dyi = —ai{t)xi — hi{t)x\ — x\ + (i^ci(t), 

where 

a^{t)=a{t) + 0{d^'^), 



h{t) = 3x^it) + 0{d), 



(90) 



(91) 
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There is a time = 0{(p/^) such that di(t) < for t < and di(t) > for t > t^,. 
Using ( |87|) and some algebra, we obtain the existence of a constant M > such that 



Ml 

M 



< < 



fd|t|-i 



M 



^ < 



fd|t|-4 

^-5/3 



for < 



(t^ -d2/3 
t > d^/\ 



(92) 



Consider the Lyapunov function 



V 



We have ^(-1) ^ Mod"^ and 



1/2 



dV 



1 ri 



21/ 



-ddiXi - hx^yi + d ciyi 



^2-'/^d\^\ + 2'/^^V\ 



From (p2|) we get the bound 



2^1/2^3 |ci(s) Ids ^ Mid^ 



Let M2 = Mo + 2Mi and define the time 



i = sup {t I y(s) < M2d'^ for -1 ^ s < t}. 

-1 ^ t sS 



(93) 



(94) 



(95) 



(96) 



For t ^ t we get from ( p^ ) and a standard result on differential inequalities (see for 
instance |Ha]| ) that 

V{t) ^ (Mo + Mi)d^ + MsMid^'^/^ < Mad^ (97) 
for sufficiently small d, which proves that t = t^, and thus F(t^) = ©(d^). 

Using the fact that di > we obtain 



dV ^ -d^V + 2-'/^d'\c,\ + 2'/'^V\ 
2 ai ai 

We choose M4 such that Mad^/^ < M4 and define 

i = sup {t I F(s) < M4d^/3 for ^ s < t}. 



For t ^t, we have 



< t < 1 



y ^ ^[l + MM4d2/3]y + 2-V2rf3|^^| 



(98) 



(99) 



(100) 



Using Gronwall's Lemma, we obtain that if MM4d2/3 ln[ai(l)/ai(t*)] ^ In 2, 

V{t) ^2\^y^\{U) + [2ai(t)]^/' fd'JSli^ds ^ M,d'l\ (101) 
.«i(i*)J A* ai(s)V2 

Taking M4 > M5, we obtain i = 1 and thus 1/(1) = 0{d^/^). Transforming back to the 
original variables, we obtain x{l) = x^{l) + 0{d^/^) and y(l) = di:^(l) + 0{d^/^). □ 
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5.4 After the Bifurcation 

In this section, we analyse equation ( |52[ ) for r ^ (e/d)^/^. It is in fact more convenient to 
use the form 

ef = T] 

(102) 

efl = fL{T)^ + 2a{T)7] + j{T)e + 6{T)^ri-f-fr, + 0{\\z\\^)+eR{^,ri,T,e), 

obtained by the transformation rj <^ rj — air)^, where fi = ^ — a?. The right-hand side 
vanishes approximately on three curves (^,r/) = (0,0) and (C,?7) = [Q-{t),Q)^ where 

ei(r) = ±/iV2 + o{t) = ±t1/2 ^ (103) 

We assume here that -R(0, 0, 0, 0) > so that the initial condition is close to , the other 
case is obtained by symmetry. The translation ^ = C+(7") + Cii ^ = yields 



(104) 



em = -Pil + 2a??i + lil + kim - - Chi + 0{\\zif) + sR, 
where (see equation ([5^)) 

a(T) = a - + i(7 - = (c - i)r + ©(r^/^), 
/?(T) = 2ef -<5e; = 2r + 0(r3/2), 

7(r) =7-2C; = -2t1/2 + 0(t), (105) 
5(r) = 5-3C; = -3x^/2 ^O(r), 

«^i(r) = -d,e;(r) = -ir-'^' + ^(l)- 

Notice that the linearization of ( |104| ) admits the eigenvalues a it i/3^/^. It can be diago- 
nalised approximately by the transformation 

C = ^ [ei-/4 /3I/4 _ e- i -/4 ^^-1/4] + 1 i ^4 ^-1/4^^ _ (log) 

V 2 V 2 

We obtain a system similar to (^): 

eC = p(t)C + e(/^(T)C + et(;(r) + (7(C, C, r, e), (107) 



where 



p(t) =a + i/3^/2^iied/3-^/^ 

V.(T) = ii/3r'-^ar'/' = 0(T-^), 

«i(r) = e-'^/^ p-^/^R{0, 0, r, 0) + ©(r^^), 



(108) 



5(C, C, r, e) = 0(r-V4|^|2) ^ o(r-i|C|=^) + ©(r-^/^KI^). 
Finally, by Corollary ^.91 , we get the following estimate on the initial condition: 

|C((e/d)2/^)| ^ (109) 
The main result of this subsection is the following. 
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Proposition 5.10. Assume that c = da/d/i(0) < ^. Then there exist positive constants 
M , n and T such that, if e and d are sufficiently small, any solution of ( |107] ) with initial 
condition ( |109| ) satisfies the hounds 

|C(r)| ^ M[er-3/4 +ei/2g-«rVe] (e/d)V3 ^ r ^ T. (110) 

Corollary 5.11. On the same time interval, there is a constant Mi > such that the 
solution of ( |102| ) satisfies 

The proof of Proposition |5.10| follows directly from the two lemmas given below. 
Lemma 5.12. Assume that c < There are positive constants M and T such that ( p!0 



admits a particular solution Co{'t) satisfying 

|Co(t)| ^ for {e/df^ ^ r s: T. (112) 



Proof: The proof is similar to the proof of Lemma |5.6| , so we only outline the differ- 
ences. 

1. We have to eliminate quadratic terms from the equation as well. In order to get a 
normal form similar to (|7lD, we start by eliminating quadratic terms, then we remove 
the term linear in C, and then only the nonresonant cubic terms. 

2. For (e/d)"^/^ ^ t ^ ^^^'^ , we eliminate the real part of the linear term by the transfor- 
mation C, = exp[cr^/e]^i. We then change the direction of time, fix a C,i{e^/'^) of order 
g5/8 g^j^^ Lgjiinia |5^. 



3. For T ^ e^/^, we proceed exactly as in Lemma 5^. □ 



Lemma 5.13. Any solution of ( |107D with initial condition ( |109D satisfies 

\C{r) - Co(r)| ^ Me^/^e— for (e/(i)i/3 ^ r ^ T. (113) 

Proof: 

• Step 1: Hamiltonian system. 

Consider, as a special case of ( |104| ) , the Hamiltonian system 

(114) 



eii = Vi- kr"'/' 



2" 

em = -2rei - SV^fi - fi- 



Lemma ^.13 shows the existence of a particular solution ^o(''") = 0{eT r/o(T) 



0{eT ^/2). If (6,m) = (Co,%) + iS,2,m), the dynamics of (6,%) is governed by a 
Hamiltonian of the form 

^(6,r/2,r) = ^4 + rh{T)e2 + V^k2{T)e^ + (115) 



where A;i(r) and fc2(T) are bounded functions. Using ( [L09D , one can show that for 
sufficiently small d, there exists a constant Mq such that 

H{t) ^ ModV^. (116) 
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Step 2: Normal forms. 
We write ( |107D in the form 



eC = p{t)C + SipirK + ewir) + g'{C, C, r) + (C, C, r, e), (117) 



where g^iCjC^T) is the contribution of the Hamiltonian approximation (114) of ( |104D , 
and 

gHC,C,r,E)=0{T'/^\C\'') + Oir-'/'\C\') + Oir-^/^\C\^). (118) 



This relation holds because the full system (|104| ) is a perturbation of size 1 + ^/T of 
the Hamiltonian system ( |114| ). 

We now perform a number of changes of variables: a translation C, = Co(''") + Cii 
where Co(''") is the particular solution of Lemma 5.12; a linear change of variables 
zi = ^2 + '5(r)^2 5 where s satisfies (|67|), which cancels the term linear in (^i, and a 
transformation to normal form C2 = Cs + ^(Cs) Csj which yields the equation 

^Cs = P3(t)C3 + c°(t; e)|C3pC3 + c\r; e)\C3\% + 0{r-^/^\C3\^), (119) 

where the functions = 0{t^^) and = 0(r^^/^) denote contributions of g^ and 
5^, respectively. 

• Step 3: Bounds on the coefficients. 
We claim that 

Rep3(r) = a(r) + ©(er-^/^)^ Rec°(T) = ©(r^^/^)^ (^20) 

The first claim can be checked by a direct calculation. We observe that the linearization 
of ( |104| ) around the particular solution (^0;??o) has the form {^p2a)^ where a = 
a + ©(er"^/^) and (5 = P + 0{eT~^^'^). Then we show that the function s(r) occuring 
in the linear transformation is such that Ims(r) = 0{eT~^). 

The second claim can be proved without lengthy calculations. By construction, tc^{t) 
is polynomial in {r^/^, er~^/^}. If we assume by contradiction that the leading term 
of Rec''(r) is of order t~^, we reach the conclusion that would grow faster than 
allowed by the estimate (116). 

• Step 4-' Final estimate. 

The Lyapunov function V = |z3p satisfies the equation 

eV - 2ktV [1 - M2r-3/V - Msr-^/V^/^j ^ (^21) 

where k = ^ — c and M2 > 0. We obtain the conclusion in a similar way as in 
Proposition □ 



6 Qualitative properties and robustness 

We conclude by discussing a few cases to which Theorem ^[2| does not apply. 

• Other values of da/dp: 

The result requires a sufficiently low value of da/d/i, which should be such that we 
traverse the bifurcation diagram Fig. ^ from region A to one of the regions D, E or 
F. Note that in a neighbourhood of the bifurcation point, the qualitative behaviour is 
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Figure 4. Same as Fig. but for values of da/d/i to which Theorem 2.2 does not apply, 
(a) When /i(r) = 2r, we move just along the C-D line of Fig. |^; the equilibria are marginally 
stable, thus the solution keeps oscillating around them with a constant amplitude, (b) 
When [i(t) = 0.5r, we go from region A to region C; now the asymmetric equilibria are 



unstable; the trajectory oscillates for some time around them because of Proposition 5.8 
but ultimately escapes to the periodic orbit. 



independent of the functions 7 and 5. The only relevant fact is that the bifurcating 
equilibrium branches should be attracting. This fact is not surprising when we steer 
into regions E or F. It is a bit more surprising when we go into region D, because the 
equilibrium is surrounded by an unstable periodic orbit. The fact that the trajectory 
lands inside this orbit is due to the specific nature of the Hamiltonian approximation 
valid near the bifurcation point. 

If we further increase da/d// so as to reach region C, the equilibria are no longer 
attracting. On the C-D boundary, trajectories oscillate around these equilibria with 
approximately constant amplitude (Fig. ^), while in the region C, they depart from 
them after a time of order e^l"^ and finally reach the outer periodic orbit (Fig. Qb). In 
view of our control problem, this behaviour is not desirable in the sense that we have 
an appearance of large amplitude oscillations, although the delay is not macroscopic. 
Imperfect control: 

While constructing the feedback control in Section ^, we had to adjust precisely three 
parameters: the imaginary eigenvalue shift 1/ and the two parameters vi and V2 of 
the nonlinear part. If these parameters are not set exactly to the desired value, the 
dynamics will still be governed by equation ([Tl|), but the functions a(T), /u(t), j{t) 
and 6{t) will not vanish exactly at the same time. The same is true to some extent if 
the feedback u{x, r) does not vanish exactly on the nominal equilibrium branch x*{t). 
In other words, such imperfections result in the fact that we traverse the bifurcation 
diagram of Fig. || on a curve which misses the origin. Two situations may occur: 
If we traverse the diagram on a path A-F-E, we first experience a pitchfork bifurcation, 
which results in an exchange of stabilities as shown in [LS|. Thus the trajectory will 
still track a stable equilibrium, but the bifurcation will occur a little bit earlier. 
If we traverse the diagram on a path A-B-C, the situation is less favourable. The Hopf 
bifurcation A-B will not be felt immediately, but when the region C is reached, there 
is a risk that the trajectory jumps to the periodic orbit. This behaviour can only be 
avoided if the slope d;^/da is large enough that the region D is reached before the 
trajectory has departed from the asymmetric branch. 
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